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ABSTRACT 

We use the Spitzer Space Telescope Enhanced Imaging Products (SEIP) and the Spitzer Archival 
Far-InfraRed Extragalactic Survey (SAFIRES) to study the spectral energy distributions of spectro- 
scopically confirmed type 1 quasars selected from the Sloan Digital Sky Survey (SDSS). By combining 
the Spitzer and SDSS data with the 2-Micron All Sky Survey (2MASS) we are able to construct a 
statistically robust rest-frame 0.1 - 100 /im type 1 quasar template. We find the quasar population 
is well-described by a single power-law SED at wavelengths less than 20 /zm, in good agreement with 
previous work. However, at longer wavelengths we find a significant excess in infrared luminosity 
above an extrapolated power-law, along with significant object-to-object dispersion in the SED. The 
mean excess reaches a maximum of 0.8 dex at rest-frame wavelengths near 100 //m. 
Subject headings: galaxies: active; quasars: general; infrared: galaxies; surveys 



1. INTRODUCTION 

Quasars, or Quasi-Stellar Objects (QSOs), are thought 
to be super-massive bla ck holes (M > 10^ Mq ) actively 
accreting material (e.g. lElvis et al]|1994|) . The energy 
output by these objec ts may control the rate of star for- 
mation in galaxies (H opkins et al.|[2006[ ). so understand- 
ing and characterizing these objects is key to galaxy evo- 
lution studies. These luminous objects are characterized 
by broad emission lines in the rest-frame optical and ul- 
traviolet and are believed to represent the peak accretion 
period for most black holes. During this period black- 
hole accretion is expe cted to dominate th e bolometric 
output of the galaxy ([Hopkins et al.|[2OO60 . and so the 
SED template should represent the energy output of an 
actively accreting black hole. 

The Spectral Energy Distribution (SED) of quasars is 
particularly important because it traces how much en- 
ergy is output by the black hole at each wavelength, and 
is needed to separate the galaxy and quasar light when 
determining key properties such as the sta r formation 
rate or stellar mass of a galaxy ^Fritz et al.l [2b06). Un- 
fortunately, spectral energy distributions in the mid- and 
far-infrared (FIR) wavelengths are poorly understood. 
To date, infrared quasar SEDs have been limited mainly 
by small samples and survey biases. 

In the absence of detailed spectrophotometry, SEDs 
can be generated by averaging broad-band pho- 
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Shang et all 120 lit iWu et al.ii201l[ ). In the mid-infrared 



(3 - 100 /x m) the most compreh ensive study to date is 
the work of iPoUetta et al.l ()2007[ ). which used an X-ray 
selected sample of 136 Active Galactic Nuclei (AGNs) 
detected at 24 /xm by the Spitzer M ultiband Imaging 
Photometer (MIPS iRieke et al.l 120041 ). However, this 
work adopted a power-law extrapolation beyond 3 /xm 
due to a lack of robust mid-infrared data. Later work by 
iHatziminaoglou et al.l (|2008.) attempted to extend these 
SEDs beyond 20 /xm using t he Spitzer Wide-are a In- 
fraRed Extragalactic f SWIRE . Eonsdale et al1l2003[ ) sur- 
vey, but reached inconclusive results due to the small 
number of high signal-to-noise detections at 70 and 160 
lim. More recent works (jHa tzimin aoglou et al. |201( ; 
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have extended quasar SEDs out to observed wavelengths 
of 500 /xm, using small samples of Herschel data. 

In this paper, we will focus on the 20 - 100 fim 
SEDs of type 1 Quasars. The goal of this work is 
to characterize this template SED and the intrinsic 
scatter around it which may indicate residual galaxy 
light and/or obscuration of the black hole. To probe 
the 3 - 100 fim range we use the Spitzer E nhanced 
Imaging Products (SEIP, ICapak et al.l lin prep l and the 
Spitzer Archival FIR Extragalactic Survey (SAFIRES, 
iHanish et aLllin prepi ) which provide science quality mo- 
saics and photo metry for most C ryogenic Infrared Array 
Camera (IRAC iFazio et al.ll2004D and MIPS data in the 
Spitzer Heritage Archive. Th e Sloan Digital Sky Survey 
(SDSS. lAbazaiian et al.|[2009l) Data Release 7 provides a 
robust SED within the observed 0.3 - 1 ^tm range, and 
the T wo Micron All-Sky Survey (2MASS, Skrutskic eraLl 
|2006[) . is used to fill in the 1-3 fim wavelength range. 

The paper is organized as follows: Section [2] describes 
the data used by the SAFIRES template creation pro- 
cess, Section 13] explains the methods used to generate our 
SED template, Section [J] presents the templates created 
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through this process, and Section [S] presents our conclu- 
sions. We assume a ACDM cosmology (fio = 0.3, 51a = 
0.7) with a Hubble constant of Hq = 70 km s~^ Mpc^^. 

2. DATA 

2.1. Parent Sample 

The process used to create an accurate quasar SED 
requires photometric data across a wide range of wave- 
lengths. We combine 15 wavelength bands ranging from 
the visible bands of the Sloan Digital Sky Survey to the 
far-infrared MIPS bands of the Spitzer Space Telescope. 
For inclusion in this work, sources were required to meet 
five criteria: 

• Inclusion in the SDSS Quasar Catalog, Data Re- 
lease 7. The optical limits associated with this cat- 
alog restricts this sample to type 1 quasars, as type 
2 quasars lack the optical luminosity to be identi- 
fied jDy the SDSS^ This is examined in more detail 
in lReves'et al.l (i2008f ) 

• A location in one of two large, well-studied ex- 
tragalactic (|b| > 20.0 degrees) fields with deep 
Spitzer coverage: the Spitzer Dee p Wide-Field 
Survey (S DWFS, 'Ashbv et al." '20091 or the Lock- 
man Hole (iLockma n ct al. 1986). 

• A SNR > 10 detection exceeding the minimum 
depth of coverage in at least one of the seven 
Spitzer bands. Only point sources were consid- 
ered; sources showing extended emission were re- 
moved from the sample. 

• A redshift of z < 3.0, provided by the SDSS Quasar 
Catalog. 

• An optical/near-infrared power-law SED slope of 
a > —1.2, removing a small number of severely 
reddened type 1 QSOs from the sample. 

Wh ile the full SDSS Quasar Catalog (ISchneider et al.l 
12010( 1 includes 105,873 objects, the first three criteria 
result in a sample of 328 type 1 QSOs possessing mea- 
surements in optical as well as mid-infrared wavelengths. 
We further refine this sample of 328 sources with two 
additional criteria, to remove small numbers of outly- 
ing sources which would have interfered with our abil- 
ity to reliably generate a characteristic type 1 template. 
High-redshift sources tend to suffer from severe blend- 
ing issues in the higher MIPS bands, and so seventeen of 
these 328 quasars are removed by our z < 3.0 criterion. 
Ten additional quasars were removed from the template 
generation algorithms by the fifth and final criterion, as 
they appeared t o be highly dust-redd ened type 1 QSOs 
as described bv iGlikman et al.l ()2012D . While these ten 
quasars are still classified as type 1 QSOs, their SEDs 
are not typical of the general type 1 population. As a 
result of our last two criteria, a total of 27 type 1 QSOs 
are removed from our SED template generation. 

These five criteria result in a sample of 301 optically 
identified quasars with measurements in at least one 
Spitzer infrared band. As shown in Figure [TJ these 
quasars show a fairly even distribution of redshifts be- 
tween 0.0 and 3.0, resulting in measurements across a 



smooth distribution of rest-frame wavelengths extend- 
ing from the Lyman Break to beyond 100 /xm in the 
far- infrared (FIR) . This consistency is important for our 
template generation process, as a more restrictive red- 
shift cutoff would result in gaps in the distribution of 
rest-frame wavelengths. As the distribution of redshifts 
in this sample of 301 quasars matches closely that of the 
SDSS sample as a whole, we feel confident that our selec- 
tion methods do not inherently bias towards or against 
quasars at certain redshifts. 

2.2. Optical and Near- Infrared: SDSS and 2MASS 

To be included in this work, a quasar requires de- 
tection in the Sloan Digital Sky Survey, Data Release 
7 (SDSS, Abazaiian et al. 2009). The SDSS Quasar 
Catalog ( Schneider et al.i ^2010) identifies 105,873 dis- 
tinct quasars within the full SDSS area, covering 9380 
square degrees of sky. Each potential quasar is con- 
firmed through spectral examination by that work; no 
explicit emission line equivalent width limit was applied 
to rule out other types of objects, and only a conserva- 
tive emission lin e FWHM limit (> 100 km s^^) was ap- 
plied. However. ISchneider et al.l ( |2010l ) use catalogs pro- 
vided by previous SDSS wo rks to explicitly re move type 
2 quasars ([Reves et al.'2008\ Seyfert galaxies (!Hao et all 
,2005.), and BL Lac objects (Plot kin et al. 2010l ). The re- 
sult is a sample believed to consist only of type 1 quasars, 
an assumption we do not dispute. 

Each SDSS quasar possesses measured luminosities in 
aU five optical SDSS wavelength bands (u', g', r', i', 
z'), with central wavelengths ranging from 350 to 850 
nm. The sample ranges from < z < 6.5 and all 
quasars are brighter than m^' < 19.1, which effectively 
removes Type 2 QSOs from our potential sample. As 
the SDSS Quasar Catalog only includes spectroscopically 
confirmed quasars, is unknown how many additional type 
1 quasars could also fall into the 'SDSS- Spitzer coverage 
region. The sky-derived fiux density limits within op- 
tical bands are consistently several orders of magnitude 
below the typical values for type 1 quasars, and so none 
of these quasars would lack SDSS data due to insuffi- 
cient flux density. The completeness of the SDSS sample 
is controlled entirely by the ability to identify quasars' 
spectra. 

Each SDSS quasar is also compared t o data from the 
Two Micron All-Sky Survey (2MASS, iSkrutskie et~all 
l2006[) . using both the 2MASS All-Sky Point Source Cat- 
alog (PSC) and the 2MASS Survey Point Source Reject 
Table. While this survey covers the entire sky, its sensi- 
tivity is far lower than that of the SDSS, preventing us 
from acquiring near-infrared fluxes for most of the SDSS 
quasars. As the available 2MASS catalogs list only point 
sources, no extended objects could be included in this 
sample; this criterion greatly reduces the potential for 
source confusion, and more closely matches the process 
used for our mid- and far-infrared bands. 

2.3. Mid-Infrared: Spitzer Enhanced Imaging Products 

(SEIP) 

The init ial release of the Spi tzer Enhanced Imaging 
Products (iCapak et al.l lin prepl ) improves upon exist- 
ing datasets in two signiflcant ways: the use of larger, 
deeper Spitzer fields, by combining contiguous data from 
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Figure 1. The redshifts for the objects included in the SDSS and SEIP/SAFIRES samples. The full SDSS catalog is shown in black, 
while the two lower histograms display all quasars located within the area observed by SAFIRES (dark grey), and those with coverage in 
either 70 fj.m or 160 ^im (light grey). Redshift bin sizes are 0.01 for the SDSS sample and 0.1 for SEIP/SAFIRES. 



multiple observation requests (AORs), and use of addi- 
tional indirect detection mechanisms to supplement the 
SNR > 10 direct measurements. This mid-infrared 
archival data is now available, and contains sources in 
the four Spitzer IRAC bands (at 3.6, 4.5, 5.8, and 8.0 
^m) as well as the 24 ^m band of MIPS. The full SEIP 
will contain sources from approximately 1500 square de- 
grees of sky, with the subset released at the present time 
spanning around 50% of the final amount. This subset 
contains over 20 million detected sources of all types, 
with an estimated 20 - 30 million additional Galactic 
sources eventually to be added. The initial release of the 
SEIP includes nine large, deep extragalactic fields; for 
this work, we select quasars from two of the largest con- 
tiguous fields, the Lockman Hole (|Lockman et al.iri986D 
and the field of the Spitzer Deep Wide-Field Survey (SD- 
WFS, Ashby ct al. 2009). 

Objects are included in the SEIP Source List if they 
were directly detected at > lOcr in any bands, and de- 
tections in multiple bands are treated as a single object 
if the detections in multiple bands have no more than 1 
arcsecond separation between IRAC bands, or 3 arcsec- 
onds when comparing to the MIPS 24 //m band. While 
301 SDSS-located type 1 quasars fall within the full area 
observed by Spitzer, only 222 - 229 quasars were located 
within the area of adequate coverage depth within each 
IRAC band. While their observed areas were equiva- 
lent, each IRAC band had a slightly different total due 
to depth of coverage differences or the effects of bad pix- 
els. 

2.4. Far-Infrared: SAFIRES 

We refine our type 1 QSO template in far-infrared 
wavelengths using source catalogs from the SAFIRES 
project. SAFIRES expands the SEIP to include mea- 
surements of point sources observed within the MIPS 70 
and 160 /xm bands throughout the cryogenic phase of 
the Spitzer mission. Unlike SEIP, SAFIRES is limited 
to only include extragalactic (|b| > 20.0 degrees) fields. 
With this limitation, and the lower total area covered 
by the MIPS 70 and 160 /im bands, the full SAFIRES 
data set will span approximately 600 deg^ of sky when 
completed. 

The methods used to extract these FIR sources are ex- 
tremely similar to those used for the MIPS 24 fim mea- 
surements, as the software used to process the cryogenic 



Spitzer data was constructed with these wavelengths in- 
cluded. The primary differences between these two MIPS 
bands and the previous five Spitzer bands involve use of 
the Germanium Reprocessing Tools (GeRT) to further 
refine data, along with an increase in the correlation dis- 
tances to compensate for the larger PSFs and pixel scales 
of the two longer-wavelength MIPS bands. Similar source 
identification mechanisms to those of the SEIP were used 
throughout the SAFIRES project. 

2.5. Detection Methods 

As explained in Section \TJ\ our core sample consisted 
of 301 type 1 quasars possessing Spitzer observations in 
at least one infrared band. Of these, 253 possess data 
within at least one MIPS band, necessary for the genera- 
tion of our far-infrared template. The number of sources 
located within the area observed by each band with suf- 
ficient depth of coverage, and the number of detections 
within those areas, are shown in Table [T] The less re- 
solved bands nearly always overlap areas observed in 
more resolved bands; less than 10% of the sources ob- 
served in MIPS 70/im or 160 /im lack coverage in MIPS 
24 /im, and almost all of the 24 /im sources possess IRAC 
coverage as well. 

We attempted to measure the magnitude of each source 
within the sample considered here, for each of the fifteen 
wavelength bands used in this work. For each measure- 
ment, there were four possible outcomes: 

• "Direct" detections. These are sources identified by 
reliable, self-contained source detection algorithms, 
using no information from other bands. For all 
Spitzer bands, a signal-to-noise ratio of at least 
10 is required, while the 2MASS catalogs required 
only a 7a detection (5cr if detected within all three 
2MASS bands). The tools used to extract this in- 
formation vary with instrumen t, with IRAC bands 
using Source Extractor (Berli n fc Arnouts I [T996[) 
and MIPS bands using the APEX package within 
MOPEX (Makovoz & Marlcau 2005). The direct 
SDSS and 2MASS values are drawn from three pro- 
vided source catalogs: the SDSS Quasar Catalog, 
the 2MASS All-Sky Point Source Catalog (PSC), 
and the 2MASS Survey Point Source Reject Ta- 
ble. Objects in the Point Source Reject Table are 
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Table 1 

SAFIRES data distributions 

Number of I SDSS 2MASS IRAC MIPS 

quasars (all bands) (all bands) 3.6/^m 4.5/^m 5.8^m S.Ofim 24fim TOfim 160/xm 

Covered 301 + 301 219 221 215 221 224 164 172 

Direct 301 * 59 * 169 173 162 171 151 11 

Indirect N/A N/A 46 43 32 36 70 89 95 

Nondetections N/A N/A 4 5 21 14 3 64 77 

Note. — f : Includes any SDSS-identified quasars within the areas observed by the Spitzer 
Space Telescope within each band. No attempt is made to quantify the number of quasars within 
this area lacking SDSS detections. 

p. All objects in the provided catalogs possessed magnitudes for all bands on an instrument; 
these were treated as direct detections. 



only used if they possessed confidence ratings above 
90%, in order to ensure reliable SED generation. 

"Indirect" detections. These consist of Spitzer 
aperture flux measurements centered on positions 
provided by the SDSS Quasar Catalog, instead of 
using positions measured directly in the band in 
question. As this process is positionally less accu- 
rate than the direct detections, larger apertures are 
used to ensure complete measurement of objects. 
This results in an increase in quoted noise levels, 
and so we require only SNR > 3. While the SDSS 
Quasar Catalog uses a similar method to add SNR 
> 2 2MASS measurements to their sources, we use 
no analogous system for our own 2MASS measure- 
ments, and there exists no similar method for the 
SDSS catalog. 

Nondetections. As the SDSS Quasar Catalog 
formed the basis for the sample considered here, 
all quasars are detected in SDSS. There are only 
a handful of true nondetections in any of the four 
IRAC bands, and only two quasars possess ade- 
quate 24 MIPS coverage depth but lack suc- 
cessful detections. While the 2MASS catalogs lack 
measurements for most of the SDSS sources, the 
majority of sources in are nondetections in these 
bands due to the prohibitive magnitude limits of 
the 2MASS projects; these magnitude limits will 
be discussed in more detail in Section [3l Nonde- 
tections form a significant part of this sample in 
the MIPS bands, for which our provided data ta- 
ble gives 3(7 upper limits, although we make more 
realistic assumptions when generating SED tem- 
plates. 

Rejections. Sources can be rejected from this sam- 
ple for a number of technical or statistical reasons, 
even if there exist image data for those regions. 
These reasons include extended features (common 
in the MIPS bands due to blending issues with 
nearby sources), cosmic rays, edge effects, and bad 
pixels. The most common reason for rejection was 
lack of coverage, as the spatial coverage of the two 
far-IR MIPS bands were smaller than those of the 
IRAC bands, where less than half of the Lock- 
man Hole region was observed in at least one of 
the two far-infrared bands. A less common issue is 
that of insufficient coverage, with individual mea- 
surements for an object ignored if the coverage in 



that band's composite image fell below a nonzero 
band-specific threshold. While most common in 
the Spitzer MIPS bands, this last effect is most 
easily seen in the IRAC bands; all IRAC-covered 
objects possessed data for all four bands, but vari- 
ations in depth of coverage occasionally led to the 
rejection of some IRAC bands for a given object, 
resulting in the differing numbers of IRAC-covered 
objects in Table [T] 

As shown in Table [U the relative contributions of these 
categories vary significantly with wavelength band. In- 
direct aperture measurements formed only a small part 
of this sample in the IRAC bands, and were in the mi- 
nority for MIPS 24 fim. Conversely, very few 70 or 160 
/im sources had direct detections, with nearly all of the 
detections in these bands using the indirect method. 

In total, 193 out of 301 sources within the SDWFS and 
Lockman Hole fields met the minimum coverage levels of 
our work in at least one FIR band, and these sources 
are essential for the construction of our far-infrared tem- 
plate. While the areas covered by the MIPS 70 and 160 
/im bands are the same, the depth of coverage minima re- 
duce the numbers of sources in both bands, with 167 pos- 
sessing adequate 70 /im coverage and 176 with 160 /im, 
as many far-infrared fields lacked the depth necessary 
for reliable measurements. The majority of these indi- 
rect far-infrared detections had directly detected counter- 
parts in at least one of the five Spitzer wavelength bands 
from SEIP. In total, only 26 of the 70 fiui and 31 of the 
160 ^m aperture detections had adequate depth of cov- 
erage in at least one mid-infrared SEIP bands while not 
possessing any direct detections, and another 21 quasars 
had sufficient depth in one of the two far-infrared bands 
but lacked coverage in MIPS 24 /im or any of the IRAC 
bands. 

Direct measurements were correlated purely by dis- 
tance, with each source matched to the nearest qualifying 
source from other catalogs. SDSS, 2MASS, and IRAC 
bands required a separation of no more than 1 arcsec- 
ond, while MIPS 24 /im detections could vary by up to 
3 arcseconds from their counterparts. The 70 /im MIPS 
sources were considered to match those of the best possi- 
ble SDSS or SEIP counterpart within a match radius of 
6 arcseconds, while the 160 /im data required a separa- 
tion of no more than 8 arcseconds. These distances were 
selected by examining the correlations for several hun- 
dred bright point source objects within an initial sample 
of quasars, and selecting the typical distance at which 
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no apparent mismatches occurred. As these correlations 
suffer greatly from confusion issues in the higher MIPS 
bands, this distance does not scale linearly with the PSF 
radius of each band. The conservative nature of this 
distance cutoff produces a correlated source list of high 
confidence, instead of a comprehensive source catalog. 

However, while the correlation distances are chosen 
such that mismatches were minimized, the large size of 
the point spread function in the 70 and 160 /xm bands 
means that the potential for source confusion is much 
highter than in a more resolved mid-infrared band, as 
any aperture is likely to include flux contributions from 
a number of nearby objects as well. The result is the dis- 
tinct possibility of a single SAFIRES detection or aper- 
ture including flux contributions from multiple distinct 
sources, potentially increasing our estimated SED at far- 
infrared wavelengths. Figure [2] illustrates the difficulties 
inherent in correlating far-infrared sources to more re- 
solved bands. As most 70 and 160 /im "sources" are ac- 
tually the combination of several distinct objects, their 
center positions are unlikely to fall within our correlation 
distances, shown as the black circle in each image. This 
effect will be examined in more detail in Section |4l 

3. SPECTRAL ENERGY DISTRIBUTIONS 

By matching the Spitzer-ohserved sources from the 
SEIP and SAFIRES projects with their counterparts in 
the SDSS/2MASS quasar catalog, as detailed in Sec- 
tion m we obtain a wide variety of data spanning up 
to fifteen bands with wavelengths ranging from 3500A to 
almost 200 /xm, with a robust population of rest-frame- 
adjusted effective wavelengths extending from the near- 
ultraviolet to around 100 fim. The measurements of each 
SAFIRES quasar for each available wavelength band arc 
given in our online table, with a format described in Ta- 
bled 

For the Lockman Hole, the areas containing sufficient 
depth of coverage in MIPS 70 and 160 fim are less than 
half of those for IRAC, resulting in a somewhat lower 
number of potential quasars from those bands. However, 
in the areas where seven Spitzer bands were available, 
the detection rates for all bands remained comparable. 
While IRAC and MIPS 24 /im allowed for a large num- 
ber of direct, SNR > 10 detections, our indirect aperture 
method was sufficient to provide SNR > 3 counterparts 
for these objects in the MIPS 70 and 160 /im bands. 
Likewise, the SDWFS field had excellent detection rates 
as the spatial coverages in all bands as the bands' respec- 
tive areas overlapped substantially more, with over 90% 
of the area in each region observed in all seven Spitzer 
wavelength bands. Within the 12.82 square degrees con- 
tained within SDWFS, the SDSS Quasar Catalog hsts 
81 quasars, of which 71 were directly detected in IRAC 
and the 24 fj,m MIPS band, all 71 of which successfully 
acquired flux density values through our indirect method 
in both the 70 and 160 ^m MIPS bands. 

The type 1 QSOs comprising this sample varied in lu- 
minosity by several orders of magnitude. As most of 
these samples had incomplete data, with several bands 
lacking adequate depth or being removed for various 
technical reasons, it was necessary to normalize each ob- 
ject's flux densities to a common baseline before any sort 
of logic could combine the data points into a coherent 
template. As we are trying to examine the shape of the 



far-infrared SED, we normalize each object based on its 
luminosities at rest-frame wavelengths between 0.2 and 
1.0 fiui, the most well-studied part of the quasar SEDs. 
Nearly every quasar (285 out of 301) possessed at least 
two measurements within this region, as every quasar in 
our sample had a complete set of five SDSS measure- 
ments, with many objects adding 2MASS detections as 
well. More than half of our sample (166 objects) possess 
at least four observations within that rest-frame wave- 
length range, resulting in a robust normalization pro- 
cess. All fluxes corresponding to rest-frame wavelengths 
within those limits were compared to the Polletta QSOl 
template value at the appropriate wavelength, to derive 
a band-specific deviation from the QSOl template; each 
object's fluxes were then normalized by the mean devi- 
ation for all bands within that rest-frame window. This 
process effectively normalizes each object's SED by its 
optical and near-infrared intensities without penalizing 
objects lacking depth in individual bands, which is es- 
sential given the wide variation in spatial coverage for 
the various bands. 

We used these normalized luminosities to generate a 
single SED templatej_with a method ba sed on that of 
several earlier works (IRic hards ct al. 2006^ IPolletta et aD 
l2007t iHatziminaoglou et al.. 2QQ8^ . Measurements are 
sorted into a series of rest-frame wavelength bins, which 
begin with a width of 0.2 dex. This size was chosen 
to minimize the chance of an object possessing multiple 
rest-frame measurements within a single bin, a very real 
possibility when dealing with the optical SDSS bands. 
We then use an iterative adjustment process to narrow 
the bins down further while keeping a minimum number 
of 100 rest-frame source measurements within each bin, 
with the final bins being as small as 0.03 dex. The re- 
sulting distribution of binned data provides greater detail 
to the templates in a wavelength regime that previously 
used a simple power-law relationship, while keeping the 
number of detections within each bin at a reasonable 
value. To minimize the impact of erroneous data points, 
we use the weighted medians within our wavelength bins 
to generate a final SED template, shown in yellow in 
Figure Oj. 

4. RESULTS 

4.1. SED Template 

We have created an SED template using the meth- 
ods described in Section El for the SAFIRES-SDSS sam- 
ple of 301 type 1 QSOs, 167 of which possess data in 
either the 70 or 160 /im wavelength bands. The re- 
sults are shown in Figure [3jj. This template, using the 
binned median of detections, is shown in ye how, while the 
previ ous SED curves ar e shown in blue (IPolletta et al.] 
[20071 QSOl) and green ([Richards et all 120061 ). We have 
also confirmed the reliability of this template through 
the use of a Kaplan-Meier ma ximum-likelihood estima- 
tor (jFeigelson fc Nelson I I1985D to account for the far- 
infrared nondetections. The Kaplan-Meier method al- 
lows for a robust estimation of the contributions of cen- 
sored data, such as nondetections caused by our band- 
specific brightness limits, but the template remains dom- 
inated by the brightest sources and so the exact method 
used to estimate luminosities for our nondetections will 
have little impact on the resulting template. Note that 
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Figure 2. Several quasars in tlie SDWFS field, shown in the IRAC 3.6 fim, MIPS 24 fim, MIPS 70 ^lm, and MIPS 160 ^lm bands. All 
images are 1 arcminute square, centered on the SDSS position for each quasar. The overlay circles show the position of each quasar, and 
the maximum correlation distances allowed within each band. 



Table 2 

SAFIRES photometry format 



Column 



Name 



Descriptions 



(1) 

(2) 

(3) 

(4) 

(5) 
(6H10) 
(11)-(13) 
(14)-(17) 
{18)-(20) 
(21)-(23) 
(24)-(27) 
{28)-(30) 



ID 
RA 
DEC 
Z 

REGID 

U.MAG, G.MAG, R.MAG, I.MAG, Z.MAG 
J.MAG, H.MAG, K.MAG 
I1.MAG, I2.MAG, I3.MAG, I4.MAG 
M1.MAG, M231AG, M3.MAG 
J.STAT, H.STAT, K.STAT 
I1.STAT, I2.STAT, I3.STAT, I4.STAT 
M1.STAT, M2.STAT, M3.STAT 



Name from SDSS Quasar Catalog 
Right Ascension (J2000) 
DecUnation (J2000) 
Spectroscopic redshift from SDSS 
Region ID from SEIP; see Capak et al. (innre^ 
Magnitudes (AB) for quasars in the SDSS u', g', r^v^a^d z' bands 
Magnitudes (AB) for quasars in the 2MASS J, H, and Kg bands 
Magnitudes (AB) for quasars in the Spitzer IRAC 3.6, 4.5, 5.6, and 8.0 fim bands 
Magnitudes (AB) for quasars in the Spitzer MIPS 24, 70 and 160 ^m bands 
Detection status for 2MASS bands: Direct or Nondetection 
Detection status for IRAC bands: Direct, Indirect, Nondectection, or Low Coverage 
Detection status for MIPS bands: Direct, Indirect, Nondetection, or Low Coverage 



Note. — The data table is available online for the 301 type 1 quasars used by the template generation algorithm. For nondetections possessing 
adequate coverage, the quoted magnitudes are the 3cr fiux density limits at the appropriate location. 



in Figure [3]i, all object and SED luminosities are nor- 
malized based on their optical and near-infrared char- 
acteristics, as described in Section [3] This counteracts 
the large inherent variation in our source measurements, 
with measured luminosity densities varying by up to five 
orders of magnitude across this sample; as the derived 
template's shape did not vary significantly between sub- 
samples divided by luminosity or redshift, we feel justi- 
fied in applying this normalization to our entire sample 
before calculating the median template. 

The shape of the SAFIRES SED template matches well 
to the SDSS, 2MASS, IRAC, and MIPS 24 ^m data used 
by other studies, confirming the shape of earlier tem- 
plates out to rest-frame wavelengths of approximately 10 
/im. The width of the filters limits the effects of individ- 
ual spectral lines; some earlier works avoided this issue 
by adding the spectral lines for a typical AGN to a binned 
median template, but we do not attempt to duplicate this 
process. While the numbers of sources identified in the 
far-infrared MIPS bands are somewhat less than their 
counterparts in the lower-wavelength IRAC and SDSS 
bands, they are sufficient to add detail to a section of 
the standard SED templates that previous studies sim- 
ply could not include. 

As our quasar sample was drawn from the SDSS 
Quasar Catalog, with its inherent optical selection bi- 
ases, no Type 2 QSOs, Seyfert galaxies, or BL Lac ob- 
jects were included within any of these samples. As ex- 



plained in Section [2.1[ ten additional dust-reddened type 
1 QSOs were removed for having optical power-law fits 
far from typical of the remainder of the sample. As a re- 
sult, the template generated by SAFIRES can be seen as 
representative of the majority of type 1 QSOs. As seen 
in FigurelSji, the SAFIRES rest-frame template generally 
matches the power-law relationship inherent to the previ- 
ous templates within the mid-infrared regime. However, 
we observe a significant deviation from these template at 
rest-frame wavelengths above 20 - 30 /xm. The template 
generated here appears to begin to deviate from the pre- 
vious QSOl template at approximately 20 fim, increasing 
to an infrared excess of approximately 0.7 dex relative to 
the QSOl templat e at rest-frame wavele ngths of around 
100 /.tm. As the IRichards et al.l ()2006l ) template pre- 
dicts a n even lower far-infrared SED than IPolletta et al.l 
(|2007f) , our results deviate from it at a wavelength closer 
to 10 /xm, and the luminosity discrepancy near 100 fim 
is more than one order of magnitude. If the Richards 
template was truly representative of type 1 quasars, few 
of our quasars would have possessed far-infrared lumi- 
nosities above the SAFIRES brightness limits. As this 
turned out to not be the case in practice, we exclude the 
Richards template from further analysis and compare our 
results only to the PoUetta et al. (2007) QSOl template. 

We wish to quantify the effect of the observed far- 
infrared excess on the total infrared luminosity of our 
quasar sample. While the classical Lm extends from 
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Figure 3. Spectral Energy Distributions of tlie SAFIRES quasar sample. 

(a) : The rest-frame SED for the sample considered here, showing the normalized, binn ed median of detection s (yellow, thick line), as well 
as normalized versions of the IPoUetta et al.i 120071') QSOl template (blue dashed) and 'Richards et al.' f2009) template (green). Colors of 
individual data points indicate the wavelength bands in which they were observed, with the SDSS bands in black, 2MASS bands in gray, 
and the seven Spitzer bands using distinct colors (violet through red), with the two SAFIRES MIPS bands denoted with square symbols. 
Arrows denote 3(t flux limits for sources with no detections despite adequate depth of coverage. 

(b) : The discrepancy between each non-normalized measurement and the QSOl template at the corresponding rest-frame wavelength. Blue 
lines denote the PoUetta TQSOl and BQSOl quartile templates. 



8 fim to 1000 nm, the SAFIRES-generated SEDs only 
extend to approximately 150 fim. Over a more lim- 
ited wavelength range, 8 - 100 yum, this template's in- 
frared flux density exceeds that of its Polletta counter- 
part by a factor of 1.53, with nearly all of the discrep- 
ancy falling above a rest-frame wavelength of 80 fim. De- 
termining the cumulative effect on IR luminosity out to 
1000 fim would require assumptions about the shape of 
the SED beyond the range of our observations. If the 
quasar SED returns to the power-law relationship of ear- 
lier works above 200 /im, then our observed discrepancy 
above 30 /im would have a negligible effect on the total 



infrared value. If this discrepancy is an indicator of a sys- 
tematic underestimation of far-IR flux in earlier quasar 
templates, however, then the total disparity could easily 
dwarf the amount quoted above. 

A similar far-infrared disparity has been observed in 
other work s extending into the far-infrared, most recently 
iDai et al.l (j2012), which used Herschel SPIRE observa- 
tions at 250, 350, and 500 /im of a small sample of 32 
quasars to estimate the mean quasar far-infrared SED. 
In that study, the mean quasar deviates from the ear- 
lier power-law templates at rest-frame wavelengths above 
~10 /im, by an average of 1.4 dex at a rest-frame wave- 
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length of 90 /im. While our observed discrepancy is 
smaller, the disparity is most likely explained by sta- 
tistical differences between th e samples and m ethodol- 
ogy of the studies in question. iDai et aP ()2012| ) removes 
many fainter or possibly blended quasars, and divides 
its sample into "FIR detected" and "FIR undetected" 
subsamples, with the quoted 1.4 dex disparity being the 
average of the former. For SAFIRES, most objects that 
would have fallen into the latter category or that would 
have been removed due to confusion issues can instead 
take advantage of our indirect methods to extract us- 
able fluxes in the absence of a direct SNR > 10 detec- 
tion. As a result, we find the two studies to be in good 
agreement, confirming a significant far-infrared excess in 
quasars when compared to the templates of earlier stud- 
ies. 

To examine these discrepancies in more detail, we com- 
pare the non-normalized quasar magnitudes to the QSOl 
template in Figure [5f), to determine the point at which 
our median SAFIRES type 1 template significantly de- 
viates from the power laws used by earlier works. The 
earlier Polletta results include not only a mean type 1 
QSO (QSOl) template, but "top quartile" (TQSO) and 
"bottom quartile" (BQSO) templates for the type 1 pop- 
ulation as well. While these quartiles are not true stan- 
dard deviations, they can provide a similar function for 
our analysis. Subtracting the QSOl template from the 
SAFIRES data, without any normalization done on in- 
dividual sources, we confirm that sources with data in 
these MIPS bands show a significant increase relative to 
the baseline template at far-IR rest-frame wavelengths, 
with medians falling well outside the curves correspond- 
ing to the QSOl quartiles in the far-infrared regime. 
Other than the 2MASS bands (shown in gray), which 
clearly suffer from their 5cr limit on point sources, each 
band shows a comparable distribution of magnitudes up 
to a rest-frame wavelength of at least 20 /im, but the 
far-infrared bands have considerably larger distributions 
of magnitudes. The mean of this distribution does not 
significantly deviate from the quartile curves until a rest- 
frame wavelength of approximately 30 /im, and the me- 
dian discrepancy reaches its maximum value of approxi- 
mately 2.0 magnitudes (0.8 dex) at around 80 /im. Few 
of our observed far-infrared data points (including only 
10 observations in the 160 /im MIPS band) fall below 
the previous Polletta QSOl template, although approx- 
imately 40% of quasars are not detected in each MIPS 
band. These nondetections are not sufficient to counter- 
act the excess in far-infrared quasar fiux demonstrated 
by our observed quasars, and so our mean far-infrared 
luminosity appears to still fall well above the QSOl tem- 
plate line. 

Using Monte Carlo logic, we have also confirmed that 
the excess relative to the QSOl template is not simply 
the statistical result of nondetections caused by the rel- 
atively high limiting fluxes in the far-infrared bands. As 
a test, we assumed all of our quasar's SEDs had fol- 
lowed the QSOl template's power-law relationship in the 
far-infrared, relative to their observed visible and near- 
infrared fiuxes. Under this assumption, more than half of 
our sources' far-infrared fiux densities should have fallen 
below the 3(t detection threshold. However, the median 
of the detected sources would average only 0.1 - 0.2 mag- 
nitudes higher than the QSOl template, far less than the 



disparity we have observed. This Monte Carlo simulation 
allows us to conclude that our observed increase in far- 
infrared flux density cannot be explained through purely 
statistical means. 

This sample contains several inherent biases for red- 
shifts, with high-redshift objects being signiflcantly less 
likely to be successfully detected in the 70 or 160 /im 
MIPS bands. The completeness of the sample within 
each wavelength band is illustrated in Table [1] while the 
majority of quasars possessed data in the IRAC band, 
the longer- wavelength MIPS bands had direct detections 
in a lower number of cases, generally losing the fainter 
high-redshift quasars within shallow Spitzer fields. This 
has the potential to raise our median template within the 
transition regions between bands, but this bias is difficult 
to quantify and does not appear to have a significant ef- 
fect on our final templates. No correction to magnitudes 
based on redshift was applied to the sample, as correctly 
applying the "K correction" requires assumptions about 
the shape of the spectrum of the object to which it is 
applied. As our goal is to establish a template for the 
shape of that spectrum, these sort of assumptions are not 
feasible. As previously mentioned, we have limited our 
SAFIRES template generation sample to include only 
those objects of redshift z < 3.0, even though the in- 
put SDSS Quasar Catalog contained QSOs with redshifts 
of up to 6.0, as shown in Figure [TJ While the regions 
possessing SAFIRES data were observed during several 
surveys of high-redshift or dusty quasars, the total pop- 
ulation distribution remains comparable to that of the 
full SDSS sample. We observed no significant changes 
in this sample when generating templates from only ob- 
jects within a specific redshift range, such as only using 
objects at z < 1.0. 

4.2. Confusion Bias 

The most significant concern when processing far- 
infrared data is source confusion. Unlike the other poten- 
tial reasons for rejection, such as bad pixels or insufficient 
depth of coverage, this has the potential to create a sig- 
nificant and systematic discrepancy in the shape of the 
SED curve. The primary cause of source confusion is tied 
to the sheer size of the Point Spread Function (PSF) in 
our less resolved bands. For MIPS 160 fj,m, the PSF has 
a width nearly seven times that of the MIPS 24 /im band, 
affecting a region approximately 40 arcseconds in radius, 
while the MIPS 70 /im PSF has a radius of around 17.5 
arcseconds. As a result, foreground objects that remain 
distinct in the more resolved bands, such as IRAC and 
the MIPS 24 /xm bands, may fall well inside the 70 /im 
and 160 /im apertures, and so be included in the mea- 
sured flux for that quasar. 

To determine whether this effect is sufficient to explain 
the noted ffux discrepancy in the far-infrared, we create 
a number of false sources in existing images, using the 
shape of the PSF for each band, at random positions 
across a single 0.25 deg^ area. Object positions were 
only constrained to stay within the area possessing our 
minimum depths of coverage, approximately half of the 
total area; no constraints were used to limit proximity 
to true sources in each image. These false objects were 
scaled to a wide range of intensities, and were designed to 
have randomly selected PSF signal-to-noise ratios rang- 
ing from 1.0 to 20.0, relative to the sky uncertainty of 
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Figure 4. The measured flux excess for an assortment of artificially created sources, relative to the intended signal-to-noise ratio of the 
false source. The dashed lines denote the 10- and 3-(t cutoffs used for direct detections and PSF/aperture measurements, respectively. The 
white lines denote binned medians within SNR bins of width 1.0 for each band. In none of the MIPS bands are the disparities between 
input and output flux large enough to explain our observed 0.8-dex discrepancy in the far-infrared SED. 

the entire iraage. We repeat this process a number of 
times to ensure statistical confidence, creating a total of 
20000 sources in MIPS 24 /.tmwithin our single image. Of 
this number, a total of 10461 false sources were located 
in areas possessing adequate depth of coverage, with the 
remaining 9539 sources removed from our analysis. Sim- 
ilar methods were used in the other MIPS bands, with 
4558 false sources placed within the same region in 70 
/im, and 2579 in the 160 /zm band. The result is shown 
in Figure 21 comparing the "true" signal-to-noise ratios 
of our created false sources to the effective increases in 
flux caused by confusion. 

The horizontal axis in Figure |4] denotes the intensity of 
the false source relative to the image's background varia- 
tion. The quoted SNR for each object is directly propor- 
tional to its total flux, but any confusion would depend 
on the intensities of the potentially confused sources near 
the created object, and so correlate more closely with a 
local flux ratio, such as SNR, than with an absolute flux. 
The vertical axis is simply the measured output flux rel- 
ative to the known intensity of the created object. High- 
SNR objects should, by definition, converge to a ratio of 
1.0, while low-SNR objects would be expected to have 
higher ratios as nearby sources contribute proportion- 
ately more to the measured flux. While there are several 
significant outliers, such as the small high- SNR popula- 
tion below the normal population in the 160 /Ltm band, 
most of these sources will not affect our actual results. 
These outliers, were they real, would have been rejected 
from the SAFIRES sample for practical reasons, such as 
a PSF extending far outside the area of spatial coverage 
within a given band. Additionally, the source measure- 
ment algorithm used for this check allows the measure- 
ment aperture to move towards local flux maxima within 
a limited radius, which would effectively increase the ef- 
fects of source confusion in cases where a nearby object is 
brighter than the source to be measured. As our quoted 
aperture fluxes do not allow for this movement, and in- 
stead are fixed at the position given by the SDSS Quasar 
Catalog, the bias quantities derived from this method 
will act as a maximum deviation. 

While this method shows a notable excess when the 
faintest sources are created, the median does not signif- 
icantly deviate from a 1.0 ratio for the vast majority of 
SAFIRES objects. The median SAFIRES SNR for all de- 



tected 160 fj,m sources is 8.54, despite every source in this 
band requiring only a Scr indirect measurement, instead 
of a lOcr direct detection, and at that signal-to-noise ra- 
tio the median increase in fiux due to confusion was less 
than 1%. Even for our faintest sources, those near the 
hard SNR > 3 limit required for fiuxes added through 
our aperture system, the median ratio does not exceed 
1.20 for 70 fim data and 1.33 for 160 /im. As a result, 
we feel confident in our assumption that the noted SED 
discrepancy cannot be explained by any sort of system- 
atic overmeasurement due to source confusion within the 
MIPS 70 and 160 fim bands. 

5. CONCLUSIONS 

We have derived a type 1 QSO SED template using 
the algorithms described in Section [Sj resulting in a ro- 
bust template capable of modeling type 1 quasar spec- 
tra across the majority of infrared, optical, and ultravi- 
olet wavelengths. The quasar SED template discussed 
here confirms the fi ndings of many recent works (e.g. 
iHatziminaoglou et al. 2008; Dai et al"1l2012| ). in that the 
spectral energy distributions of quasars do not follow a 
single power law beyond wavelen gths of 3 lim, as wa s 
assumed in earlier works such as iPolletta et al.l ()2007f) . 
Assuming that the sample used in this work is not biased 
towards any specific type of type 1 quasar, the SAFIRES 
template indicates a significantly higher flux above a rest- 
frame wavelength of 20 /im when compared to previ- 
ous templates, potentially resulting in a signiflcantly in- 
creased total far-infrared luminosity. Within this sample, 
it appears that while all type 1 quasars have similar SEDs 
in the visible and near-infrared regimes (barring the ten 
removed from our initial sample) , their far- infrared char- 
a cteristics can vary dras tic ally. 

IHatziminaoglou et al.l (|2010( ) and iBonfield et al.l 
()2011|) explained the possible excesses in their far- 
infrared MIPS measurements through a combination 
of stellar emission and a cold dust component linked 
to bursts of star formation in an extended accretion 
disk. Similar r esults were seen in the recent work of 
iDai et aP (|2012l ). where the projected dust temperatures 
needed to explain the observed far-infrared excess were 
typical associated with starburst galaxies. Whether 
our observed discrepancy could be better modeled by 
further separating the type 1 QSO category, defined 
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primarily by AGN characteristics, into raultiple distinct 
subtypes with differing disk star formation properties 
is beyond the scope of this work. It is possible that 
our observations could indicate a number of subclasses 
within the greater type 1 quasar category, with the 
various subclass' differences only becoming apparent at 
infrared wavelengths of more than 100 /xm as the other 
characteristics of the host galaxies begin to dominate. 

Other far-infrared telescopes co uld be used to improve 
this work; many recent works (Hatziminaoglou et aLl 
2010: Bonfield et al. 2011; Dai et al. 2012; Sajina et aT| 
2OI2I) use data from the Herschel SPIRE instrument to 
examine quasars further into the far-infrared, allowing 
them to better characterize the downturn in quasar SEDs 
above a rest-frame wavelength of 100 /im. The Herschel 
PACS instrument covers wavelengths similar to those of 
MIPS, albeit with greatly improved sensitivity and res- 
olution, and so can perform an analysis similar to that 
shown in this work. However, the primary constraints 
on the analysis in this work are those tied to the sam- 
ple selection criteria in Section [2?T1 such as the available 
far-IR coverage area and inclusion in the SDSS quasar 
catalog. The spatial resolution and far-infrared depth of 
coverage of the instruments are less important than the 
total number of quasars observed. 

As the Spitzer Enhanced Imaging Products and its 
SAFIRES offshoot continue to process and release 
archival data through the Spitzer Heritage Archive, the 
number of quasars with corresponding archival Spitzer 
data will rapidly increase. The SEDs shown in this work 
will be updated with progressively larger numbers of 
sources as new data permits. Additionally, the SDSS 
Quasar Catalog used in this project was based on SDSS 
Data Release 7; as the corresponding list for SDSS DR9 
is expected to increase the number of verified quasars 
by almost 70%, we would expect a comparable increase 
in the number of quasars included in our analysis. 
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